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Abstract - We show gapped critical environment could remarkably prevent an enhanced decay 
of decoherence factor and quantum correlations at the critical point, which is nontrivially different 
from the ones in a gapless critical environment (Quan, et.al Phys. Rev. Lett. 96, 140604 (2006)). 
The quantum correlations display very fast decaying to their local minimum at the critical point 
while maximum decaying occurs away from this point. In particular, our results imply that collapse 
of decoherence factor is not indicator of a quantum phase transition of environment as opposed to 
what happens in a gapless criticallity. In the week coupling regime, the relaxation time, at which 
the quantum correlations touch rapidly local minima, shows a power-law singularity as a function 
of gap. Furthermore, quantum correlations decay exponentially with second power of relaxation 
time. Our results are important for a better understanding and characterisation of gap critical 
environment and its ability as entanglers in open quantum systems. 
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Quantum correlations (QCs) are of primary importance 
in quantum information and quantum computation 

They are related to the basic issue of understand¬ 
ing the nature of non-locality in quantum mechanics [6][^ . 
Quantum systems used in quantum information process¬ 
ing inevitably interact with the surrounding environment. 
These correlated surrounding systems induce quantum de¬ 
coherency which plays a key role in the understanding of 
the quantum to the classical transition [^[^. As a re¬ 
sult, in the last decade a lot of efforts have been devoted 
to investigate QCs dynamics and decoherence factors of 
central systems in various environments 10 11 . The de¬ 


the short time region the LE decays as a Gaussian. How¬ 
ever for long time limits they found that it approaches an 
asymptotic value, which strongly depends on the strength 
of the transverse magnetic held. Further, the decoherence 
of a system coupled to a spin environment with QPT has 
been investigated 
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coherence of a system coupled to a spin environment with 
quantum phase transition (QPT) has been investigated in¬ 
tensively in various studies [Tl}j20]. Quan et al. 21 


con- 


The quantum phase transition occurs at a level cross¬ 
ing point (gapless phase transition) or converged point of 
avoided level crossing point (gapped critical point) [^ . 
Because of the convergence of the energy levels at the criti¬ 
cal point (CP), some special dynamic features may appear 
in the dynamic evolution of the central system in contact 
with an environment with QPT. It is shown that the dis¬ 
entanglement of the central system is greatly enhanced by 
the gapless quantum criticality of the environment 12 


sidered induction of the Ising-type correlated environment 
on the Loschmidt echo (LE), and found that the decaying 
behavior of LE is best enhanced by gapless QPT of the 

surrounding system. Rossini et al. 22 depicted that in there exist separate sta tes^ sho wing non-classical behaviour 


Furthermore, the decoherence induced by the critical en¬ 
vironment may display some universal features 13 . Since 


without entanglement 24-26 , the quantum entanglement 


P-1 


















R. Jafari and A. Akbari 


Entanglement: A=3 Discord: A=3 Entanglement: A=0 



0 4 8 12 16 20 0 4 8 12 16 20 0 4 8 12 16 20 

t t t 


Fig. 1: (Color online) (a) Density plot of entanglement as a function of time, for Jo = 1, and Je = 4. (b) The variation of 
the quantum discord as function of time t where the Hamiltonian parameters set as Jo = 1, and Je = 4. (c) Entanglement 
evaluation for the gapless critical environment (Jo = Je = 1). We consider N = 400 and weak-coupling between qubits and 
spin bath ^ = 0.1. 


does not include all type of quantum correlations and 
their non-classical properties. Therefore, other measures 
of quantum correlations are expected. 


A new but promising notion is quantum discord first 
introduced by Ollivier and Zurek 
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which can effec¬ 
tively capture all QCs in a quantum system. Recently, 
the dynamics of quantum discord under the effect of the 
environmental spin chain has increasingly been investi¬ 
gated 28 . The results show that quantum discord is 


more robust than entanglement for the system subjected 
to the environment 


29 


In this letter, we consider two spins coupled to the one¬ 
dimensional general quantum compass model (GQCM) 
in the presence of a transverse field [34}|^. 


GQCM is a simplified model that describes the nature of 


the orbital states in the case of a twofold degeneracy 34 


We analyze the effect of the gapped critical environment 
on the dynamic evolution of the two-spin entanglement, 
quantum discord, and negativity. Our comprehensive 
results show that induction of the gapped critical envi¬ 
ronment is nontrivially different than the gapless critical 
environment 13 21 , and strongly recommend it as a 
better entangler and quantum channel. In particular, 
in the weak coupling regime, QCs at the gapped CP 
decay rapidly from a maximum and after a very fast 
initial transient, start oscillating around an average value 
whereas show maximum decay away from CP. The results 
show that at CP an average value of QCs is enhanced 
with the increasing of the gap. 


The full Hamiltonian which considers two spins cou¬ 
pled transversely to the spin 1/2 general quantum compass 


chain is fully characterized by 77 = + 77/ with 12 39 


N' 

'^E ^ + Je^2i ^ ^2i+l + + '^2i) 


i=l 


N' 




( 1 ) 


i=l 


with 77/ describing the interaction between the central 
two qubits and the surrounding spin chain with 

coupling strength g. Additionally, He is the Hamilto¬ 
nian of the environment describing the one-dimensional 
GQCM (Id-GQCM). In this representation, Id-GQCM 
is constructed by antiferromagnetic order of X and Y 
pseudo-spin components on odd and even bonds at which 
the pseudo-spin operators are constructed as linear combi¬ 
nations of the Pauli matrices (cr“='^’*'’^): = ai{±0) = 

cos Oaf ± sinffcrf 38 . Here 9 (—9) is arbitrary angle rel¬ 
ative to cr® for even (odd) bounds. Jg and Jo charac¬ 
terise the even and odd bound couplings respectively, h is 
a transverse field, and N = 27V' is the number of spins. 

We should emphasise that the Id-GQCM is exactly solv¬ 


able with the Jordan-Wigner transformation 34,38,40 
which in momentum space leads to 


He = 




( 2 ) 


where ( 7 ^’'^) denote independent fermions creation 
(annihilation) operators. For states with even fermions 

77/ = Va -I- \/b and 77/ = \/a — '/h, with a = 47i^ -|- | Jp -I- 
|L|2 and b = {16h^ + 2 |J|2)|L|2 + where 

the parameters L and J are defined by L = (Jq-I-J ee*^)/4, 
and J = (Joe'®-Jee‘('=-®))/4. 

It is known that the ground state is separated from the 
lowest-energy pseudo-spin excitation by a pseudo-spin gap 
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First we con- 


which vanishes at cos 0 c = h/y/JoJe 
centrate on an idiosyncratic case of 0 c = 7''/2 in the ab¬ 
sence of external magnetic field. For this case, the ground- 
state has a macroscopic degeneracy of away from 

the isotropic point {Jo ^ Je), which becomes 2 ^^^ when 
the orbital interactions are isotropic. The gap at fc = tt 
is given by the anisotropy of the pseudo-spin exchange, 
A = |Je — Jo|, and just vanishes at Jg = Jo which im¬ 
plies that the degeneracy increases by an additional factor 
of 2 due to the band-edge points. It has been proven 
that in the absence of an external magnetic field, for 
0c = 7r/2, GQCM with Z 2 symmetry is critical for ar¬ 
bitrary Je/Jo 39 . QPT takes place between two dif¬ 


ferent disordered phases where the model exhibits highest 
possible frustration of interactions 38 39 . The isotropic 


point corresponds to a multicritical point where the gap 
closes quadratically at fc = ±7r as a result of the confluence 


of two Dirac points 41 


Since the central two qubits operators {o'%<Jb) 
the environmental spin chain {uf) satisfy the commu¬ 
tation relation [a\ + cr|j,(T“] = 0, the total Hamilto¬ 
nian can be rewritten as H = \ ® 

|(/?^), {fi = 1,..., 4) denotes the fith eigenstate of the op¬ 
erator §(< 75 . 'J^b) corresponding to the /rth eigenvalue 
(ei = g, £2 = £3 = 0, £4 = —g). Here T-L^ is defined 
through "Hb by replacing h with , where = h e^. 
We suppose that the initial state of the total system is 
disentangled with ptot{Q) = Pab{Q) ® Pe{Q)- Consid¬ 
ering I'0 b) as an initial state of the environmental spin 
chains, and pABiO) and pe{0) = |'0b)(0b| as initial den¬ 
sity matrix state of the two-qubits system and environ¬ 
ment respectively, the evolution of the total system will 
be governed by ptot{t) = U{t)ptot{0)UHt). Accordingly, 
the reduced density matrix of two-qubits AB is obtained 
by tracing out the environment [12|, 



Fig. 2: (Color online) (a) The scaling behavior of the relaxation 
time Tr in terms of energy gap for N = 400, Jo = 1, Je = 4 
and g — 0.1. (b) Evaluation of entanglement as a function of 
time for different energy gap values, (c) The scaling behaviour 
of relaxation for small values of gap. 


written as 


/! + , 


PAB{t) = 


V 


0 

1 - 

c* 

0 


1 - 
0 


C/3 

0 
0 

1 + Cz y 


(4) 


with cp = {cx - Cy)Fi 4 , and = {c^ + Cy)F 23 . 

To quantify the entanglement dynamics of two qubits 
AB, we utilize the concurrence directly to calculate the 
entanglement of formation (EoF). The entanglement is 
a monotonically increasing function of concurrence {C'), 
which is defined by 44 


4 

PAB{t) = F^y,{t){ip,,\pAB{0)\>Py)\>p^){>Py.\, (3) 


where the decoherence factors are achieved by F^^y = 
( 0 b|G^"V^ 10 b), and {t) = is the 

time evolution operator driven by . By considering 
the ground state of the environment spin chain as an ini¬ 
tial state, the decoherence factors reduce to a LE form 
given in Ref. 21 , which is a dynamical version of the 


ground-state fidelity 40,42 


We assume that the two qubits AB initially stem from 
the A-structure states pab(O) = \{Iab + J2a Coo'd 0 cr%) 
with the identity operator on two-qubits system, Iab ■ 
The parameters Ca=x,y,z are chosen to be real that insure 
that pab{0) is a legal quantum state. This state is cho¬ 
sen in a general form to contain Bell-diagonal states and 
Werner states. According to Eq. ([^, the reduced density 
matrix in the standard basis (|tt)i Iti)) lit)) lit)) can be 


EoF = -fc log 2 fc - (1 - /c') log 2 (l - fc), 


with fc' = |(1 -I- \/l — C"2) and C = max{Ai — A 2 — 
A 3 — A 4 , 0}, where A^ are the square roots of the eigen¬ 
values in descending order of the operator pAB{t){cr\ (g) 
c’b)Ptb(^)(c’a ® o"^)- One can directly calculate the 
concurrence for the state defined by Eq. Q as C" = 
max{ ^ Moreover, quantum discord 

is given by 


Qab = |c.^|) log 2 (l - Cz ± |cT,|)-b 

(1 -b Cz ± |C; 3 |) l 0 g 2 {l + Cz ± \cp\) - c(^pAB{t)^ 


(5) 


where C{pab) = i i?), is the classical 

correlation with •§ = max{|cz|, }. To investi¬ 

gate the time evolution of QCs, we set the parameters 
Cx,z = —Cy = 1 , so that the initial state becomes the 
Bell state (| ii) -b | \\))/y/2. For this case, concur¬ 
rence and quantum discord become C = |Fi 4 (t)| and 
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Q = J2± log 2 (l ± 1 ^ 14 ( 01 )) respectively. To il¬ 

lustrate the dynamical properties of QCs, we carry out 
the numerical calculation using the exact expression. The 
density plots of the time evolution of entanglement and 
quantum discord have been depicted in the weak coupling 
limit in Fig. [^a & b). We set the Hamiltonian param¬ 
eters in such a way that the system becomes gapped at 
CP 9c = 7r/2. Obviously both entanglement and quantum 
discord are monotonically increasing functions of the de¬ 
coherence factor |Fi 4 (t)|. Thus, quantum discord behaves 
in a similar way as entanglement does and is always less 
than entanglement in the process of evolution. 

In contrast with the above situation, Fig.j^c) represents 
a density plot of entanglement for the gapless critical case. 
As is clear, the enhanced decay of QCs induced by quan¬ 
tum criticality of the surrounding environment is broken 
by gapped quantum critical environment and the maxi¬ 
mum decaying happens away from CP. The result shows 
that in the weak coupling regime, QCs decay from their 
maximum and after an initial transient, start oscillating 
around an average value. As observed in Fig. QCs are 
symmetric with respect to 9 = tt/2 and the numerical cal¬ 
culation shows that the valley narrows as g decreases and 
system size, N, increases. A more detailed analysis also 
shows that the relaxation time at which QCs decay 
to their local minimums at CP of the environment (see 
Fig. [^b)), reveals a power-law singularity as a function 
of the gap A. This is presented in Fig. [^a) which speci¬ 
fies a linear behaviour of In(Tr) versus ln(A). The scaling 
behaviour is obtained as with the exponent 

T = —0.75 for small value of the gap whereas the exponent 
T equals —0.85 for large values of the gap. Thus, one can 
conclude that the appearance of the two energy scales in 
the system is appropriate with the number of states which 
are involved. It means that the dynamics of QCs at very 
large energy gap values mainly originates from the ground 
state and the excited states have a very tiny contribution. 
It would be worth mentioning that the interaction cou¬ 
pling g does not affect the exponent r in the weak coupling 
limit, as displayed by numerical simulations. 

To study the scaling behaviour of QCs at CP, we have 
derived the scaling behaviour of local minimum values of 
QCs versus the second power of the relaxation time. This 
has been plotted in Fig. [^a), which shows the linear be¬ 
haviour of ln[EoF(Tr)], the same as ln[(5(T'r)], versus T^. 
In other words, the local minimum values of QCs scale ex¬ 
ponentially with the second power of the relaxation time, 
EoF(Tr) = exp(—with exponent S oc {—A/g). In 
particular, our results imply that the decay of the local 
minimum of QCs at CP of environment enhances with 
decreasing energy gap and increasing interaction coupling 
(see Fig. |^b)). 

However, oscillations of QCs around an average value 
increases as the energy gap of the environmental spin 
chain increases. On the other hand. Fig. [^c) reveals 
an interesting phenomenon in QCs at CP {9c = 7r/2). 



Fig. 3: (Color online) (a) Scaling of the minimum of quantum 
correlations at critical point {0^ = 7r/2) for systems of various 
gaps (Fig. (b) Entanglement evaluation as a function of 

time for different energy gap for N = 400, Jo = 1, and g = 
0.05. (c) Entanglement as a function of time at critical point 
{9c = 7r/2) for various system sizes, and Jo = 1, Je = 4, and 
5 = 0 . 1 . 


The periods of the revival of QCs is independent of the 
size of the environment which means that CP is a scale 
invariance point where quantum fluctuations extend over 
all length scales. Moreover, a similar result also has 
confirmed the above argument for the time evolution of 
QCs from mixed states Cx = 1, and Cy = —Cz € [0,1]. 
The presented results presented here indicate that in 
the strong coupling regime decoherence factors and QCs 
decay to zero in a very short time as opposed to what 
happens in the weak coupling regime. As the coupling 
strength increases, the valleys widen and the influence of 
the energy gap on the generation of QCs decreases even 
when approached along the gapped critical point. 


In the presence of a magnetic field, GQCM shows a 
gapless phase transition and our results support previous 
studies 12-^. We also examine two qutrits coupled to 
one Id-GQCM |^. The negativity displays almost the 
same dynamical behaviour as entanglement and quantum 
discord do. 


In summary, using the general quantum compass model 
as an environmental system, the dynamical evolution of 
the decoherence factors, quantum correlations, and neg¬ 
ativity of the central spins has been investigated for dif¬ 
ferent initial states. The relation between the quantum- 
classical transition of the central system, and the occur¬ 
rence of an avoided level crossing quantum phase transi¬ 
tion in its surrounding system has been analysed. It is well 
known that the gapless quantum criticality enhances de¬ 
caying of decoherence factors, while our calculations repre- 
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sent a different story for gapped critical environment. The 
finding results show that long-time quantum correlations 
at the critical point is an effect of gapped criticality, and 
maximum decaying occurs away from the critical point. 

The role of the gapped critical spin chain is to prevent 
the complete drain of information from central systems 
to the environment 46 and provides them a better 
environment for preserving quantum correlations. In 
other words, the amount of decoherence which travels 
into the central spin state depends on the excited states 
of the environment. Hence the energy gap could block the 
propagation of decoherence along the environment and 
consequently reduces its effect on the central spin. These 
results highlight the current outlook of using quantum 
spin chains as entanglers or quantum channels in quantum 
information devises 47 48 . Besides, quantum gapped 


criticality may have potential applications in quantum 
computations. 
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